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1 Introduction 

Understanding charmonium production in hadronic collisions is a challenge for quantum chro- 
modynamics (QCD). The ]/ip production cross section measurements at the Tevatron 111 |2l 
were found to disagree by about a factor of 50 with theoretical color-singlet calculations f3\. 
Soon after, the CDF experiment reported a Xci/ Xci cross section ratio that extended up to 
Pt{]/^) — lOGeV/c, where pj is the transverse momentum, and favored Xci production over 
Xc2 EI- The cross section ratio was also studied recently at the Large Hadron Colllider (LHC) 
in Ref. fS]. These measurements independently suggest that charmonium production carmot 
be explained through relatively simple models. 

This paper presents a measurement of the Xo.1 XcX cross section ratio by the Compact Muon 
Solenoid (CMS) experiment at the LHC in pp collisions at a center-of-mass energy of 7TeV. 
The measurement is based on the reconstruction of the Xc radiative decays to J/t/; + 7, with the 
low transverse momentum photons (less than 5 GeV/ c) being detected through their conversion 
into electron-positron pairs. The analysis uses data collected in 2011, corresponding to a total 
integrated luminosity of 4.6 fb^^. When estimating acceptance and efficiencies, we assume 
that the Xci and XcX are produced unpolarized, and we supply the correction factors needed to 
modify the results for several different polarization scenarios. 

Due to the extended reach in transverse momentum made possible by the LHC energies, the 
cross section ratio measurement is expected to discriminate between different predictions, such 
as those provided by the fcx"factorization |6| and next-to-leading order nonrelativistic QCD 
(NRQCD) [7J theoretical approaches. 

The strength of the ratio measurement is that most theoretical uncertainties cancel, including 
the quark masses, the value of the strong coupling constant ^s, as well as experimental uncer- 
tainties on quantities such as integrated luminosity, trigger efficiencies, and, in part, reconstruc- 
tion efficiency. Therefore, this ratio can be regarded as an important reference measurement to 
test the validity of various theoretical quarkonium production models. With this paper, we 
hope to provide further guidance for future calculations. 

2 CMS detector 

A detailed description of the CMS apparatus is given in Ref. (SI. Here we provide a short 
summary of the detectors relevant for this measurement. 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diam- 
eter. Within the field volume are the silicon pixel and strip tracker, the crystal electromag- 
netic calorimeter and the brass /scintillator hadron calorimeter. Muons are measured in gas- 
ionization detectors embedded in the steel return yoke. In addition to the barrel and endcap 
detectors, CMS has extensive forward calorimetry. 

The inner tracker measures charged particles within the pseudorapidity range \r]\ < 2.5, where 
fj = — ln[tan(0/2)], and 9 is the polar angle measured from the beam axis. It consists of 1440 
silicon pixel and 15 148 silicon strip detector modules. In the central region, modules are ar- 
ranged in 13 measurement layers. It provides an impact parameter resolution of ~15 Jim. 

Muons are measured in the pseudorapidity range \ri\ < 2.4, with detection planes made using 
three technologies: drift tubes, cathode strip chambers, and resistive plate chambers. Match- 
ing the muons to the tracks measured in the silicon tracker results in a transverse momentum 
resolution between 1 and 1.5%, for pj values up to 50 GeV/c. 
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4 Event reconstruction and selection 



The first level (LI) of the CMS trigger system, composed of custom hardware processors, uses 
information from the calorimeters and muon detectors to select the most interesting events. The 
high-level trigger (HLT) processor farm further decreases the event rate from around 100 kHz 
to around 300 Hz, before data storage. The rate of HLT triggers relevant for this analysis was 
in the range 5-10 Hz. We analyzed about 60 million such triggers. 

3 Experimental method 

We select Xci srid Xc2 candidates by searching for their radiative decays into the J/j/' + 7 final 
state, with the ]/xp decaying into two muons. The XcO has too small a branching fraction into 
this final state to perform a useful measurement, but we consider it in the modeling of the sig- 
nal lineshape. Given the small difference between the ]/tp mass, 3096.916 ± 0.011 MeV/c^, and 
the Xci and Xc2 masses, 3510.66 ± O.OZMeV/c^ and 3556.20 ± 0.09 MeV/c^, respectively [9J, the 
detector must be able to reconstruct photons of low transverse momentum. In addition, ex- 
cellent photon momentum resolution is needed to resolve the two states. In the center-of-mass 
frame of the Xc states, the photon has an energy of 390 MeV when emitted by a Xci and 430 MeV 
when emitted by a Xc2- This results in most of the photons having a pj in the laboratory frame 
smaller than 6 GeV/ c. The precision of the cross section ratio measurement depends crucially 
on the experimental photon energy resolution, which must be good enough to separate the two 
states. A very accurate measurement of the photon energy is obtained by measuring electron- 
positron pairs originating from a photon conversion in the beampipe or the inner layers of the 
silicon tracker. The superior resolution of this approach, compared to a calorimetric energy 
measurement, comes at the cost of a reduced yield due to the small probability for a conver- 
sion to occur in the innermost part of the tracker detector and, more importantly, by the small 
reconstruction efficiency for low transverse momentum tracks whose origin is displaced with 
respect to the beam axis. Nevertheless, because of the high Xc production cross section at the 
LHC, the use of conversions leads to the most precise result. 

For each Xci,2 candidate, we evaluate the mass difference Am = nififij — niiifi between the 
dimuon-plus-photon invariant mass, m^i^, and the dimuon invariant mass, nl^^. We use the 
quantity Q = Am + niy^, where niy^ is the world-average mass of the ]/ip from Ref. O, as 
a convenient variable for plotting the invariant-mass distribution. We perform an unbinned 
maximum-likelihood fit to the Q spectrum to extract the yield of prompt Xci and Xc2 as a func- 
tion of the transverse momentum of the J/t/?. A correction is applied for the differing accep- 
tances for the two states. Our results are given in terms of the prompt production ratio Rp, 
defined as 

^ ^ ^(PP ^ Xc2 + X)g(^c2 ^ ]/^ + 7) ^ ^Xc2 £i 

P (7(pp ^ ;tci + X) 5 [Xci ^y^ + l) N^^, ■ £2 ' 

where cr(pp Xc + ^) are the Xc production cross sections, B{Xc ]/^ + 7) are the Xc 
branching fractions, N^^. are the number of candidates of each type obtained from the fit, and 
ei/£2 is the ratio of the efficiencies for the two Xc states. The branching fractions B{Xci,2 — ^ 
J/i/' + 7), taken from Ref. [9], are also used to calculate the ratio of production cross sections. 

4 Event reconstruction and selection 

In order to select Xc signal events, a dimuon trigger is used to record events containing the 
decay J/t/J — ?> The LI selection requires two muons without an explicit constraint on 
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their transverse momentum. At the HLT, opposite-charge dimuons are reconstructed and the 
dimuon rapidity y{}i}i) is required to satisfy \y{}i}i) \ < 1.0, while the dimuon px rnust exceed a 
threshold that increased from 6.5 to 10 GeV/c as the trigger configuration evolved to cope with 
the instantaneous luminosity increase. Events containing dimuon candidates with invariant 
mass from 2.95 to 3.25 GeV/c^ are recorded. Our data sample consists of events where multiple 
pp interactions occur. At each bunch crossing, an average of six primary vertices is recon- 
structed, one of them related to the interaction that produces the Xc in the final state, the others 
related to softer collisions (pileup). 

In the ]/'ip selection, the muon tracks are required to pass the following criteria. They must 
have at least 11 hits in the tracker, with at least two in the pixel layers, to remove background 
from decays-in-flight. The per degree of freedom of the track fit must be less than 1.8. To 
remove background from cosmic-ray muons, the tracks must intersect a cylindrical volume 
of radius 4 cm and total length 70 cm, centered at the nominal interaction point and with its 
axis parallel to the beam line. Muon candidate tracks are required to have pj > 3.3GeV7c, 
I?/ 1 < 1.3 and match a well-reconstructed segment in at least one muon detector llTOl . Muons 
with opposite charges are paired. The two muon trajectories are fitted with a common vertex 
constraint, and events are retained if the fit probability is larger than 1%. If more than one 
muon pair is found in an event, only the pair with the largest vertex x^ probability is selected. 
For the final Xci ^rid Xc2 selection, a dimuon candidate must have an invariant mass between 
3.0 and 3.2GeV/c2 and < 1.0. 

In order to restrict the measurement to the prompt ]/ip signal component, the pseudo-proper 
decay length of the ]/tp (i]/^), defined as iy^ = Lxy ■ niy^/ pj{]/tp), where Lxy is the most prob- 
able transverse decay length in the laboratory frame IjTTl , is required to be less than 30 fim. 
In the region < 30 fim, we estimate, from the observed distribution, a contamination 
of the nonprompt component (originating from the decays of B hadrons) of about 0.7%, which 
has a negligible impact on the total systematic uncertainty. 

To reconstruct the photon from radiative decays, we use the tracker-based conversion recon- 
struction described in Refs. |IT2l - [r4J . We summarize the method here, mentioning the further re- 
quirements needed to specialize the conversion reconstruction algorithm to the Xc case. The al- 
gorithm relies on the capability of iterative tracking to efficiently reconstruct displaced and low 
transverse momentum tracks. Photon conversions are characterized by an electron-positron 
pair originating from a common vertex. The e+e^ invariant mass must be consistent with zero 
within its uncertainties and the two tracks are required to be parallel at the conversion point 
and to separate from each other only in the transverse plane. 

Opposite-sign track pairs are first required to have more than four hits and a normalized x^ 
less than 10. Then the reconstruction algorithm exploits the conversion-pair signature to dis- 
tinguish between genuine and misidentified background pairs. Information from the calorime- 
ters is not used for conversion reconstruction in our analysis. The primary pp collision vertex 
associated with the photon conversion, see below, is required to lie outside both track helices. 
Helices projected onto the transverse plane form circles; we define cfm 

as the distance between 

the centers of the two circles minus the sum of their radii. The value of dm is negative when the 
two projected trajectories are intersecting. We require the condition —0.25 < < 1.0 cm to be 
satisfied. 

In order to reduce the contribution of misidentified conversions from low-momentum dis- 
placed tracks that are artificially propagated back to the silicon tracker, the two candidate con- 
version tracks must have one of their two innermost hits in the same silicon tracker layer. 
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5 Acceptance and efficiencies 



The distance along the beam line between the extrapolation of each conversion track candidate 
and the nearest reconstructed event vertex must be less than five times its estimated uncer- 
tainty. Moreover, among the two event vertices closest to each track along the beam line, at 
least one vertex must be in common 

A reconstructed primary vertex is assigned to the reconstructed conversion by projecting the 
photon momentum onto the beamline and choosing the closest vertex along the beam direc- 
tion. If the value of the distance is larger than five times its estimated uncertainty, the photon 
candidate is rejected. 

The primary vertex associated with the conversion is required to be compatible with the recon- 
structed j/ip vertex. This requirement is fulfilled when the three-dimensional distance between 
the two vertices is compatible with zero within five standard deviations. Furthermore, a check 
is made that neither of the two muon tracks used to define the ]/xp vertex is used as one of the 
conversion track pair. 

The e+e^ track pairs surviving the selection are then fitted to a common vertex with a kinematic 
vertex fitter that constrains the tracks to be parallel at the vertex in both the transverse and 
longitudinal planes. The pair is retained if the fit probability is greater than 0.05%. If a track 
is shared among two or more reconstructed conversions, only the conversion with the larger 
vertex probability is retained. 

Only reconstructed conversions with transverse distance of the vertex from the center of the 
mean pp collision position larger than 1.5 cm are considered. This requirement suppresses 
backgrounds caused by track pairs originating from the primary event vertex that might mimic 
a conversion, such as from Dalitz decay, while retaining photon conversions occurring 
within the beampipe. 

Finally, each conversion candidate is associated with every other conversion candidate in the 
event, and with any photon reconstructed using calorimeter information. Any pairs of con- 
versions or conversion plus photon with an invariant mass between 0.11 and 0.15GeV/c^, cor- 
responding to a two-standard-deviation window around the mass, is rejected. We have 
verified that the rejection requirement, while effectively reducing the background, does not 
affect the Rp measurement within its uncertainties. 

Converted photon candidates are required to have px > 0.5GeV/c, while no requirement is 
imposed on the pseudorapidity of the photon. 

The distribution of the photon conversion radius for Xc candidates is shown in Fig.jl] The first 
peak corresponds to the beampipe and first pixel barrel layer, the second and third peaks corre- 
spond to the two outermost pixel layers, while the remaining features at radii larger than 20 cm 
are due to the four innermost silicon strip layers. The observed distribution of the conversion 
radius agrees well with MC simulations l|T4 | , which indicates a good knowledge of the material 
in the tracking volume. 

5 Acceptance and efficiencies 

In the evaluation of Rp, we must take into account the possibility that the geometric acceptance 
and the photon reconstruction efficiencies are not the same for Xci ^rid Xc2- 

In order to determine the acceptance correction, a Monte Carlo (MC) simulation sample of 
equal numbers of Xci ^rid Xc2 has been used. This sample was produced using a PYTHIA IITSi 
single-particle simulation in which a Xci or Xc2 is generated with a transverse momentum dis- 
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Figure 1: Distribution of the conversion radius for the Xc photon candidates. 

tribution produced from a parameterized fit to the CMS measured tp{2S) spectrum 1161 . The 
use of the ip{2S) spectrum is motivated by the proximity of the tp{2S) mass to the states under 
examination. The impact of this choice is discussed in Section |7| 

Both Xc states in the simulation are forced to decay to J/i/^ + 7 isotropically in their rest frame, 
i.e., assuming they are produced unpolarized. We discuss later the impact of this assumption. 
The decay products are then processed through the full CMS detector simulation, based on 
Geant4 l|T7l[T8l , and subjected to the trigger emulation and the full event reconstruction. In 
order to produce the most realistic sample of simulated Xc decays, digitized signals from MC- 
simulated inelastic pp events are mixed with those from simulated signal tracks. The number 
of inelastic events to mix with each signal event is sampled from a Poisson distribution to 
accurately reproduce the amount of pileup in the data. 

The efficiency ratio £i/e2 for different ]/xp transverse momentum bins is determined using : 

£1 

^2 NiT^NiT' 



where N^^"^ is the number of Xc candidates generated in the MC simulation within the kine- 
matic range |i/(J/t/')| < 1-0, ^1(7) > 0.5GeV/c, and N^^^ is the number of candidates recon- 
structed with the selection above. The resulting values are shown in Table [ij where the uncer- 
tainties are statistical only and determined from the MC sample assuming binomial distribu- 
tions. The increasing trend of £i/e2 is expected, because pj{]/tp) is correlated with the pj of 
the photon, and at higher photon pj our conversion reconstruction efficiency is approximately 
constant. Therefore, efficiencies for the Xci arid the Xc2 are approximately the same at high 

This technique also provides an estimate of the absolute Xc reconstruction efficiency, which 
is given by the product of the photon conversion probability, the Xc selection efficiency, and, 
most importantly, the conversion reconstruction efficiency, which corresponds to the dominant 
contribution. This product varies as a function of pi (7)/ and goes from 4 X lO^^ at 0.5 GeV/c to 
around 10^^ at 4 GeV/c, where it saturates. 
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7 Systematic uncertainties 



Table 1: Ratio of efficiencies t\ltx as a function of the J/j/? transverse momentum from MC 
simulation. The uncertainties are statistical only. 



PT(J/i/?)[GeV/c] El/ £2 



7-9 0.903 ± 0.023 

9-11 0.935 ± 0.019 

11-13 0.945 ± 0.021 

13-16 0.917 ± 0.022 

16-20 0.981 ± 0.031 

20-25 1.028 ± 0.049 



6 Signal extraction 



We extract the numbers of Xc\ ^rid Xci events, ^Xc\ ^xa' respectively, from the data by 
performing an unbinned maximum-likelihood fit to the Q spectrum in various ranges of J/i/? 
transverse momentum. 

Because of the small intrinsic width of the Xc states we are investigating, the observed signal 
shape is dominated by the experimental resolution. The signal probability density function 
(PDF) is derived from the MC simulation described in Section |5| and is modeled by the super- 
position of two double-sided Crystal Ball functions |19| for the Xc\ ^rid Xci and a single-sided 
Crystal Ball function for the XcQ- Each double-sided Crystal Ball function consists of a Gaussian 
core with exponential tails on both the high- and low-mass sides. We find this shape to pro- 
vide an accurate parameterization of the Q spectra derived from MC simulation. When fitting 
the data, we fix all the parameters of the Crystal Ball function to the values that best fit our 
MC simulation and use a maximum-likelihood approach to derive '^xa ^^'^ '^xa' which are 
the integrals of the PDFs for the two resonances. The background is modeled by a probability 
distribution function defined as 

NMcg(Q) = (Q-^or-e(e-'/o)-ft, 

where ci\ and f)\ are free parameters in the fit, and c\q is set to 3.2 GeV/c^. 

In Fig. [2] we show the Q distribution for two different ranges, 11 < piQ/i/') < 13GeV/c (left) 
and 16 < pxQ/V') < 20GeV/c (right). This procedure is repeated for several ranges in the 
transverse momentum of the J/t/^ in order to extract '^^a ^^'^ '^xa the corresponding bin. 
The results are shown in Table |2j where the reported uncertainties are statistical only. 

Table 2: Numbers of Xc\ arid Xci events extracted from the maximum-likelihood fit, and the 
ratio of the two values. Uncertainties are statistical only. 



PT(J/t/^)[GeV/c] N;,, l^xa ^xJ^x.. 

7-9 618 ±31 315 ±24 0.510 ± 0.049 

9-11 1680 ±49 788 ±37 0.469 ± 0.027 

11-13 1819 ± 51 819 ± 38 0.451 ± 0.025 

13-16 1 767 ± 51 851 ± 39 0.482 ± 0.027 

16-20 1269 ±43 487 ± 30 0.384 ± 0.028 

20-25 642 ±31 236 ± 22 0.368 ± 0.040 



7 Systematic uncertainties 

Several types of systematic uncertainties are addressed. In particular, we investigate possible 
effects that could influence the measurement of the numbers of Xc\ arid Xci from data, the 



7.1 Uncertainty from the mass fit and Xci ai^d Xc2 counting 
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Figure 2: The distribution of the variable Q = m^^^^ — m^^ + my^, for Xc candidates with pj{]/ip) 
ranges shown in the figures. The Hne shows the fit to the data. 

evaluation of £i/e2 from the MC simulation, and the derivation of the Rp ratio. In Table |3]the 
various sources of systematic uncertainties and their contributions to the total uncertainty are 
summarized. The following subsections describe how the various contributions are evaluated. 

Table 3: Relative systematic uncertainties on Rp for different ranges of ]/xp transverse momen- 
tum from different sources and the total uncertainty. 



PtO/iP) range [GeV/c ] 


7-9 


9-11 


11-13 


13-16 


16-20 


20-25 


Source of uncertainty 




Systematic uncertainty (%) 




Background shape 


1.4 


1.5 


0.9 


1.2 


1.8 


2.4 


Signal shape 


1.4 


3.0 


1.1 


1.5 


1.5 


2.3 


Simulation sample size 


2.6 


2.0 


2.2 


2.4 


3.1 


4.8 


Choice of priXc) spectrum 


4.5 


3.7 


2.9 


1.9 


0.6 


1.1 


Total uncertainty 


5.5 


5.4 


3.9 


3.6 


4.0 


5.9 



7.1 Uncertainty from the mass fit and Xci and Xci counting 

The measurement of the ratio N^^^ I ^Xci could be affected by the choice of the functional form 
used for the maximum-likelihood fit. The use of an alternative background parameterization, 
a fourth-order polynomial, results in systematically higher values of the ratio N^^^^ I ^xa' while 
keeping the overall fit quality as high as in the default procedure. From the difference in the 
numbers of signal events using the two background parameterizations, we assign the system- 
atic uncertainty from the background modeling shown in Table |3] 

We evaluate the systematic uncertainty related to the parameterization of the signal shape by 
varying the parameters derived from the MC simulation within their uncertainties. The results 
fluctuate within 1-3% in the various transverse momentum ranges. We assign the systematic 
uncertainties from this source, as shown in Table |3] 

The method to disentangle and count the Xci and Xc2 states is validated by using a PYTHIA 
MC simulation sample of inclusive ]/ip events, including those from Xc decay, produced in 
pp collisions and propagated through the full simulation of the detector. The ratio N^^^ / N^^^ 
derived from the fit to the Q distribution of the reconstructed candidates in the simulation 
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7 Systematic uncertainties 



is consistent with the actual number of Xc events contributing to the distribution, within the 
statistical uncertainty, for all ]/tp momentum ranges. Therefore, we do not assign any further 
systematic uncertainty on the determination of N;^^^ I ^xa ■ 

The stability of our analysis as a function of the number of primary vertices in the event has 
been investigated. The number of Xc candidates per unit of integrated luminosity, once trigger 
conditions are taken into account, is found to be independent of the instantaneous luminosity, 
within the statistical uncertainties. In addition, the measured ratio N^^^/N^^^ is found to be 
constant as a function of the number of primary vertices in the event, within the statistical 
uncertainties. Thus, no systematic uncertainty due to pileup is included in the final results. 

7.2 Uncertainty on the ratio of efficiencies 

The statistical uncertainty on the measurement of £i / e2 from the simulation, owing to the finite 
size of the MC sample, is taken as a systematic uncertainty, as shown in Table |3] 

Since the analysis relies on photon conversions, the effect of a possible incorrect simulation of 
the tracker detector material is estimated. Two modified material scenarios, i.e., special detector 
geometries prepared for this purpose, in which the total mass of the silicon tracker varies by up 
to 5% from the reference geometry, are used to produce new MC simulation samples |20|. With 
these models, local variations of the radiation length with respect to the reference simulation 
can be as large as +8% and —3%. No significant difference in the ratio of efficiencies is observed 
and the corresponding systematic uncertainty is taken to be negligible. 

Several choices of the generated Pt{Xc) spectrum are investigated. In particular, the use of the 
measured ]/rp spectrum IfTH gives values that are compatible with the default tp{2S) spectrum 
used for the final result. The choice of the spectrum affects the values of £i / £2 only inasmuch as 
we perform an average measurement in each bin of pT(]/xp), and the size of these bins is finite. 
We choose to assign a conservative systematic uncertainty by comparing the values of ei/ei 
obtained with the 1/^(25 ) spectrum with those obtained in the case where the pj{Xc) spectrum 
is taken to be constant in each pj bin. The corresponding systematic uncertainties are given in 
Table i 

7.3 Xc polarization 

The polarizations of the Xci arid Xc2 are unknown. Efficiencies are estimated under the assump- 
tion that the two states are unpolarized. If the Xc states are polarized, the resulting photon an- 
gular distribution and transverse momentum distributions will be affected. This can produce 
a change in the photon efficiency ratio £i/£2- 

In order to investigate the impact of different polarization scenarios on the ratio of the effi- 
ciencies, we reweight the unpolarized MC distributions to reproduce the theoretical Xc angular 
distributions ||2lJ|22l for different Xc polarizations. We measure the efficiency £i/£2 for the Xci 
being unpolarized or with helicity m^^^ = 0, ±1, in combination with the Xc2 being unpolar- 
ized or having helicity m^^^ = 0/ ^2 in both the helicity and Collins-Soper [.23] frames. The 
ratio of efficiencies for the cases involving nip^^^ = ±1 is between the cases with m^^^ = and 
m^^2 = i2. Tables |4] and |5] give the resulting £i/£2 values for each polarization scenario in dif- 
ferent ]/xp transverse momentum bins for the two frames, relative to the value of the ratio for the 
unpolarized case. These tables, therefore, provide the correction that should be applied to the 
default value of £i/£2 in each polarization scenario and each range of transverse momentum. 



7.3 Xc polarization 
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Table 4: The efficiency ratio £1/62 for different polarization scenarios in which the x^x is either 
unpolarized or has heHcity rtij^^^ = 0, ±1 and the Xo. is either unpolarized or has helicity m^^^ = 
0, ±2 in the helicity frame, relative to the unpolarized case. 



Polarization scenario {m^^-i^ r'^xa) 7 — 9 


9-11 


pT(J/i/')[GeV/c] 
11 - 13 13 - 16 


16-20 


20-25 


(Unpolarized, 0) 


0.89 


0.87 


0.85 


0.86 


0.85 


0.86 


(Unpolarized, ±2) 


1.20 


1.20 


1.21 


1.20 


1.20 


1.17 


(0, Unpolarized) 


0.83 


0.84 


0.85 


0.85 


0.85 


0.86 


(±1, Unpolarized) 


1.08 


1.07 


1.07 


1.07 


1.07 


1.07 


(0,0) 


0.74 


0.73 


0.72 


0.73 


0.72 


0.74 


(0,±2) 


1.00 


1.01 


1.03 


1.02 


1.02 


1.01 


(±1,0) 


0.95 


0.93 


0.91 


0.97 


0.90 


0.92 


(±1,±2) 


1.29 


1.29 


1.29 


1.28 


1.28 


1.25 



Table 5: The values of £i/£2 for different polarization scenarios in the CoUins-Soper frame, 
relative to the impolarized case. 



Polarization scenario {ni^^^,ni^^^ 7 — 9 


9-11 


pT(J/t/')[GeV/c] 
11-13 13-16 


16-20 


20-25 


(Unpolarized, 0) 


1.04 


1.06 


1.08 


1.07 


1.08 


1.08 


(Unpolarized, ±2) 


0.97 


0.95 


0.93 


0.93 


0.92 


0.92 


(0, Unpolarized) 


1.04 


1.05 


1.06 


1.07 


1.07 


1.06 


(±1, Unpolarized) 


0.98 


0.97 


0.97 


0.96 


0.96 


0.97 


(0,0) 


1.08 


1.12 


1.14 


1.15 


1.16 


1.14 


(0,±2) 


1.01 


0.99 


0.98 


0.99 


0.98 


0.98 


(±1,0) 


1.02 


1.03 


1.04 


1.04 


1.04 


1.04 


(±1,±2) 


0.95 


0.92 


0.90 


0.90 


0.89 


0.89 
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9 Summary 



7.4 Branching fractions 

The measurement of the xo. to Xc\ production cross section ratio is affected by the uncertainties 
on the branching fractions of the two states into J/i/? + 7. The quantity that is directly accessible 
in this analysis is Rp, the product of the ratio of the Xci to Xc\ cross sections and the ratio of the 
branching fractions. 

In order to extract the ratio of the production cross sections, we use the value of 1.76 ± 0.10 for 
-^(/tci ]/^ + l')/S{Xc2 — > ]/'^ + 7) as derived from the branching fractions and associated 
uncertainties reported in Ref. HI. 



The results of the measurement of the ratio Rp and of the ratio of the Xc2 to Xci production cross 
sections for the kinematic range Pt(7) > 0.5GeV/c and lyQ/t/')! < 1-0 are reported in Tables|6] 
and|7| respectively for different ranges of pj{]/ip). The first uncertainty is statistical, the second 
is systematic, and the third comes from the uncertainty on the branching fractions in the mea- 
surement of the cross section ratio. Separate columns are dedicated to the uncertainty derived 
from the extreme polarization scenarios in the helicity and Collins-Soper frames, by choosing 
from Tables |4] and |5] the scenarios that give the largest variations relative to the unpolarized 
case. These correspond to i^Xd'^Xci) — (^1/^2) and {mxciri^xa) — (0/0) for both the helicity 
and Collins-Soper frames. Figure |3] displays the results as a function of the J/i/? transverse mo- 
mentum for the hypothesis of unpolarized production. The error bars represent the statistical 
uncertainties and the green bands the systematic uncertainties. 

We compare our results with theoretical predictions derived from the /cj-factorization f6l and 
NRQCD [7J calculations in Fig. |4] The A:x-factorization approach predicts that both Xci and 
Xci are produced in an almost pure helicity-zero state in the helicity frame. Therefore, in our 
comparison, we apply the corresponding correction on the ratio of efficiencies from Table |4j 
amounting to a factor of 0.73, almost independent of pj. The theoretical calculation is given in 
the same kinematic range (pt(7) > 0.5GeV/c, |i/(J/'/')| < 1.0) as our measurement. There is 
no information about the Xc polarization from the NRQCD calculations, so we use the ratio of 
efficiencies estimated in the unpolarized case for our comparison. The prediction is given in the 
kinematic range Pt(7) > OGeV/c, |i/(J/i/')| < 1.0. We use the same MC simulation described 
in Section |5] to derive the small correction factor (ranging from 0.98 to 1.02 depending on pj, 
with uncertainties from 1 to 4%) needed to extrapolate the phase space of our measurement to 
the one used for the theoretical calculation. The uncertainty on the correction factor stemming 
from the assumption of the Xc transverse momentum distribution is added as a systematic 
uncertainty. The values of Rp after extrapolation are shown in Table |8] The comparison of 
our measurements with the A;T-factorization and NRQCD predictions are shown in the leftand 
rightplots of Fig. |4j respectively. The A;T-factorization prediction agrees well with the trend of 
Rp versus transverse momentum of the J/t/?, but with a global normalization that is higher by 
about a factor two with respect to our measurement. On the other hand, the NRQCD prediction 
is compatible with our data within the experimental and theoretical uncertainties. 



8 Results and discussion 
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^(PP ^ Xci + ^)B{Xci ]/^ + 7) 
^(PP ^ Xci + ^)B{Xci ]/^ + 7) 
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Table 6: Measurements of for the given pj{]/tp) ranges in the fiducial kinematic 

region pt (7) > 0.5GeV/c, \y{]/ip)\ < 1.0, assuming unpolarized ;j(;c production. The first un- 
certainty is statistical and the second is systematic. The last two columns report the additional 
uncertainties derived from the extreme polarization scenarios in the helicity (HX) and Collins- 
Soper (CS) frames. 



pj{]/xp)[GeV/c] 


0-(,Xc2)l3(Xc2) 


HX 


CS 


7-9 


0.460± 0.044 (stat.) ± 0.025 (syst.) 


+0.136 
-0.121 


+0.037 
-0.023 


9-11 


0.439± 0.025 (stat.) ± 0.024 (syst.) 


+0.128 
-0.119 


+0.052 
-0.035 


11-13 


0.426± 0.024 (stat.) ± 0.017 (syst.) 


+0.125 
-0.117 


+0.059 
-0.042 


13-16 


0.442± 0.025 (stat.) ± 0.016 (syst.) 


+0.125 
-0.121 


+0.065 
-0.044 


16-20 


0.377± 0.028 (stat.) ± 0.015 (syst.) 


+0.106 
-0.104 


+0.059 
-0.042 


20-25 


0.379 ± 0.041 (stat.) ± 0.022 (syst.) 


+0.094 
-0.097 


+0.055 
-0.040 



Table 7: Measurements of cr{xc2)/ c{Xci ) for the given pj{]/xp) ranges derived using the branch- 
ing fractions from Ref. |9|, assuming unpolarized Xc production. The first uncertainty is sta- 
tistical, the second is systematic, and the third from the branching fraction uncertainties. The 
last two columns report the uncertainties derived from the extreme polarization scenarios in 
the helicity (HX) and Collins-Soper (CS) frames. 



pjO/xp)[GeV/c] t^(;\:c2)/(^(;\:ci) hx cs 



7-9 


0.811±0.078 (stat.) ± 0.045 (syst.) ± 0.046(BR) 


+0.239 
-0.213 


+0.066 
-0.041 


9-11 


0.774±0.044 (stat.) ± 0.042 (syst.) ± 0.044(BR) 


+0.225 
-0.209 


+0.092 
-0.061 


11- 


-13 


0.752±0.042 (stat.) ± 0.029 (syst.) ± 0.043(BR) 


+0.221 
-0.207 


+0.105 
-0.074 


13- 


-16 


0.78±0.044 (stat.) ± 0.028 (syst.) ± 0.044(BR) 


+0.221 
-0.213 


+0.115 
-0.078 


16- 


-20 


0.665±0.049 (stat.) ± 0.027 (syst.) ± 0.038(BR) 


+0.187 
-0.184 


+0.104 
-0.074 


20- 


-25 


0.669±0.072 (stat.) ± 0.039 (syst.) ± 0.038(BR) 


+0.165 
-0.172 


+0.096 
-0.070 



Table 8: Measurements of ^^e given pj(J/xp) ranges after extrapolating the mea- 

surement to the kinematic region Pt(7) > and assuming unpolarized Xc production. The 
first uncertainty is statistical and the second is systematic. The last column reports the largest 
variations due changes in the assumed Xc polarizations. 



pT{]/xp)[GeV/c] 


0'(Xc2)B{Xc2) 

t^(;tci)6fei) 


Polarization 


7-9 


0.451 ± 0.043 (stat.) ± 0.025 (syst.) 


+0.137 
-0.153 


9-11 


0.427 ± 0.024 (stat.) ± 0.023 (syst.) 


+0.134 
-0.144 


11-13 


0.421 ± 0.024 (stat.) ± 0.017 (syst.) 


+0.133 
-0.142 


13-16 


0.441 ± 0.025 (stat.) ± 0.017 (syst.) 


+0.138 
-0.143 


16-20 


0.365 ± 0.027 (stat.) ± 0.016 (syst.) 


+0.114 
-0.115 


20-25 


0.387 ± 0.042 (stat.) ± 0.026 (syst.) 


+0.109 
-0.105 
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9 Summary 
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Figure 3: Ratio of the Xci to Xci production cross sections (circles) and ratio of the cross sections 
times the branching fractions to J/t/^ + 7 (squares) as a function of the ]/tp transverse momentum 
with the hypothesis of unpolarized production. The error bars correspond to the statistical 
uncertainties and the green band corresponds to the systematic uncertainties. For the cross 
section ratios, the 5.6% uncertainty from the branching fractions is not included. 
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Figure 4: Comparison of the measured values with theoretical predictions from the 

/cj-factorization ||6l (left) and NRQCD [7J (right) calculations (solid red lines). The error bars 
and green bands show the experimental statistical and systematic uncertainties, respectively. 
The measurements in the leftplot use an acceptance correction assuming zero helicity for the 
Xc, as predicted by the Aij-factorization model. The measurements in the rightplot are corrected 
to match the kinematic range used in the NRQCD calculation and assume the Xc are produced 
unpolarized. The measurements assuming two different extreme polarization scenarios are 
shown by the long-dashed blue and short-dashed green lines in the rightplot. The 1-standard- 
deviation uncertainties on the NRQCD prediction are displayed as the dotted red lines. 
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as a function of the ]/ip transverse momentum up to pj{]/ip) = 25GeV/c for the kinematic 
range Pt(7) > 0.5GeV/c and lyQ/j/')! < 1-0 in pp collisions at i/s = 7TeV with a data sample 
corresponding to an integrated luminosity of 4.6 fb^. The corresponding values for the ratio 
of the Xc2 to Xci production cross sections have been determined. 

The results have also been shown after extrapolating the photon acceptance down to zero pj. 
The effect of several different Xc polarization scenarios on the photon reconstruction efficiency 
has been investigated and taken into account in the comparison of the experimental results 
with two recent theoretical predictions. This is among the most precise measurements of the 
Xc production cross section ratio made in hadron collisions, and extends the explored ]/ip pj 
range of previous results. These measurements will provide important input to and constraints 
on future theoretical calculations of quarkoniiun production. 
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